Power law Kohn anomalies and the excitonic transition in graphene 
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Dirac electrons in graphene in the presence of Coulomb interactions of strength /? have been shown 
to display power law behavior with j3 dependent exponents in certain correlation functions, which 
we call the mass susceptibihties of the system. In this work, we first discuss how this phenomenon 
is intimately related to the excitonic insulator transition, showing the explicit relation between 
the gap equation and response function approaches to this problem. We then provide a general 
computation of these mass susceptibilities in the ladder approximation, and present an analytica l 
computation of the static exponent within a simplified kernel model, obtaining rjo = ^1 — /3//3c . 
Finally we emphasize that the behavior of these susceptibilities provides new experimental signatures 
of interactions, such as power law Kohn anomalies in the dispersion of several phonons, which could 
potentially be used as a measurement of j3. 



I. INTRODUCTION 

One of the most intriguing aspects of grapheneA contin- 
ues to be the physics of the electron-electron interaction 
and its effects on the electronic properties. This very 
active field of research^ is now of great relevance as cur- 
rent experiments are displaying unambiguous signatures 
of interactions'^''*. 

What makes this problem so interesting, among other 
reasons, is the fact that graphene, -when undoped, is a 
rather special material 'with respect to the Coulomb in- 
teraction. Because the lo'w-energy quasiparticles of this 
system are massless Dirac fermions, "with vanishing den- 
sity of states at the Fermi level, the 1/r Coulomb in- 
teraction remains long ranged even in the presence of 
screening. This makes this system very different from 
materials "with a finite Fermi surface, and gives rise to 
very unusual many-body effects. This behavior stems 
from the fact that the Hamiltonian of Coulomb interact- 
ing Dirac fermions does not contain an intrinsic energy 
scale. The kinetic and interaction terms scale in the same 
"way, and the Hamiltonian is in fact scale invariant. In 
the RG sense, the Coulomb interaction is marginal in this 
system-. The strength of the Coulomb interaction, the 
fine structure constant of graphene, may be expressed as 
a dimensionless number /3 = e^/{evp), where vp is the 
Fermi velocity of the Dirac quasiparticles, e their charge, 
and e the dielectric constant of the substrate. While a 
naive estimate gives an upper limit /3 2 for e — 1, 
implying that Coulomb interactions should be relatively 
important, experimentally their strength is still a matter 
of debate^i^. 

The presence of this interaction has many non-trivial 
consequences. From a weak coupling perspective, the 
first to be predicted was the logarithmic renormaliza- 
tion of the Fermi velocity^, which has been recently 
measured-. Other observables have been shown to be 
affected by a similar renormalizatioEk2i£. An electronic 
inverse lifetime that is linear instead of quadratic with 
energjjl is also a characteristic feature of this system. 

Another interesting phenomenon induced by interac- 
tions is the possibility of a phase transition to a broken 



symmetry state in the strong coupling regime. For a 
strong enough value of the coupling /?, low energy elec- 
trons and holes may bind into excitons, opening a gap in 
the system. The broken symmetry state is thus known 
as the excitonic insulator, and the characteristics of this 
transition have been studied throughly^^l^. The critical 
coupling is thought to be of order /? 1 but its precise 
value is still under discussiony^iSj and experimentally an 
intrinsic gap in graphene has not been observed^. This 
excitonic transition can also be considered as the many- 
body counterpart of the supercritical screening of an ex- 
ternal Coulomb potential by Dirac fermion o^^d^ . 

A more recent prediction of an interaction induced 
phenomenon, also rooted in the scale invariance of this 
system, is the presence of power law behavior with 
interaction dependent exponents in certain correlation 
functionsi^ii^"— , as if the system were in a critical phase. 
These correlators can be considered the mass-mass re- 
sponse functions (or mass susceptibilities) of the system. 
In analogy to the charge-charge response function, which 
measures the charge expectation value as a response to 
an external perturbation coupling to the charge, this sus- 
ceptibility is the analog quantity built with mass ver- 
tices (i.e. the matrices that gap the Dirac spectrum) in- 
stead of charge ones. The experimental detection of these 
correlations would represent a new signature of interac- 
tion effects beyond those already discussed, and would 
be highly desirable. In addition, the singularities of the 
mass susceptibilities are directly related to instabilities 
to spontaneous mass generation, and therefore their ob- 
servation will shed light on the problem of the excitonic 
transition. This observation is however challenging, as 
there is no simple experiment to directly probe a mass 
(as opposed to charge or spin) susceptibility. 

The objective of this work is twofold. First we will 
show the explicit correspondence between the physics of 
the excitonic transition and the behavior of the mass sus- 
ceptibilities. We will illustrate it in detail by comparing 
the transition as seen from a gap equation and a response 
function perspective. Second, we will argue that there 
are in fact current experiments available to observe the 
mass susceptibilities. We will see that there are clear 
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experimental signatures of their presence in the disper- 
sion relation of certain phonons, which acquire power law 
Kohn anomalies with /? dependent exponents. In addi- 
tion, we will show that the static exponent of the mass 
susceptibility, obtained within perturbative RG— i^Si^i, 
or numericalljii^'i^, can in fact can be obtained ana- 
lytically from the ladder approximation and is given by 

The paper is organized as follows. In Sec. |TT]we review 
the physics of the excitonic transition and its description 
in terms of a gap equation. In Sec. IIIII we describe the 
general mass-mass response function and compute it in 
different approximations. In Sec. IIVI we discuss how the 
mass susceptibility may be observed in the dispersion 
relation of particular phonons, and in Sec. |V]we present 
our conclusions. 



II. THE EXCITONIC TRANSITION 
A. The model 



The Kekule and CDW masses are time reversal invari- 
ant, and may be grouped into the three components of a 
spin 1/2 vector 



Mn — iT,j:ax,Tya,j:,az), 



(3) 



because they transform as such under the SU(2) valley 
symmetry. Therefore, when one of them is generated, the 
valley symmetry is spontaneously broken to U(l), analo- 
gously to a ferromagnet. The Haldane mass, on the other 
hand, is a scalar under valley symmetry but breaks time 
reversal. In this work we will be concerned only with 
the time reversal invariant masses and valley symmetry 
breaking. The structure of these masses and the symme- 
tries they break is analogous to the one found in QED3, 
where the valley symmetry is known as chiral symmetry. 
In fact, the problem of chiral symmetry breaking and 
mass generation in QED3 in the 1 /N approximation^ii^ 
has much in common with the excitonic transition, the 
main difference being that the interaction in QED3 is a 
Lorentz invariant gauge field (rather than an instanta- 
neous charge-charge interaction). 



We start by reviewing the excitonic insulator transi- 
tion in graphene. We will consider spinless graphene for 
simplicity, as spin will not play any role and can be ac- 
counted for when necessary. In graphene, the low energy 
excitations around the K and K' points can be modeled 
in terms of a Dirac Hamiltonian in two spatial dimensions 



H = ivp id r^\axdx + OLydy)^) 



(1) 



with d — {Tzax, cTy), where the a and r matrices act on 
the sublattice and valley degrees of freedom, respectively. 
The Coulomb interaction is included as 



Hint 2 



(2) 



This Hamiltonian has an SU(2) valley symmetry gener- 
ated by the matrices r„ = {TxCry,Tyay,Tz), in the sense 
that the SU(2) rotation ip — )• e*"^"^"-0 leaves the Hamil- 
tonian invariants^. 

When (i is large enough, this system has an instability 
towards the pairing of electrons and holes. This pairing 
generates a mass for the Dirac fermions, opening a gap in 
the spectrum, and so this state is known as the excitonic 
insulator. This instability is signaled by the development 
of a finite expectation value of a mass operator {^ip''Mipy 
There are four possible mass matrices M for the Hamilto- 
nian ([I} , which correspond to different microscopic mech- 
anisms for the instability, and which break different sym- 
metries. The first of them is the generation of a charge 
imbalance between sublattices (a charge density wave), 
and corresponds to a matrix cr^. The next two are pro- 
duced by a bond density wave, the Kekule distortior^, 
and correspond to the matrices Txcr^ and Tya^- Finally, 
a pattern of circulating currents in the unit cell gives rise 
to the mass CTzT^, known as the Haldane mass. 



B. The gap equation 

The most common approach to show the existence of 
the excitonic transition is by means of a self-consistent 
gap equation^"^"*. This is equivalent to the Hartree-Fock 
or mean field approximation. We will now review the 
main features of this approach, with the aim of high- 
lighting its relation with the response function approach 
to be discussed in the next section. For simplicity we 
set vp — I henceforth and thus the Coulomb coupling 

To derive the gap equation, the full electron propagator 
G is expressed in terms of the self-energy as = Gq^ — 
S, with the bare propagator 



kp + dk 
'^"i^) = 172 To — 



kl-k"^ + ie' 



(4) 



The self-energy is a 4x4 matrix that satisfies the 
Schwinger-Dyson equation 



-t I -^rO(fc,p)H(p - k) (Go 1 - S) \ (5) 



(27r)3 



where G, H are the full electron and photon propagators, 
and r*^ the full Coulomb vertex. This equation is so far 
exact, but to solve it one needs to assume approximate 
forms for H and F^. The simplest approximation is to 
assume a bare Coulomb vertex F'' = 1 and a bare photon 
n = 2Trf3/\q\, which gives 

This equation can be seen as the resummation of the 
"rainbow" diagrams, depicted diagrammatically in Fig. 
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[TJ To solve this equation S can be expanded in a basis 
of all sublattice and valley matrices. We will make the 
further approximation that only terms proportional to 
mass matrices are important (neglecting Fermi velocity 
and wavefunction renormalizations) i.e. we assume that 



(7) 



where summation over repeated indices is always im- 
plicit. Projecting into each mass channel (i.e. taking 
the trace of Eq. (jS]) with Mm) we arrive at 



d^k A„(fc) 



(2^)2 fc| 



fr2 



A„(fc)2 



(8) 

While the gap equations for the three masses can be stud- 
ied separately, we see that they are in fact related by 
symmetry, so that we can drop the index n. Because the 
interaction is kg independent in this approximation, the 
gap is too, and we may integrate over fco to obtain 



A(p) = TTp 



(fk A(fc) 



1 



(2^)2 fc| + A(fc)2 



(9) 



This is the simplest version of the gap equation for the 
excitonic problem. More refined approximations have 
been considered in the literature, for example at finite 
temperature^., with the renormalization of the fermion 
spectrun*i2iiiiii^, or with statioi^iii and dynamioi^ RPA 
(1/N) screening for the interaction. In general, all ap- 
proaches agree that there is an excitonic transition, but 
critical couplings vary significantly. For our discussion 
we retain the simplest form given by Eq. 



p-k 



+ 



P 



FIG. 1. Self consistent diagrammatic equation for the electron 
Green function (full line) in the rainbow approximation 



To make further progress, the standard approach is to 
implement an expansion of the gap in circular harmonics. 



^("p) — 



2^ 



'-A(p), 



(10) 



and keeping only the first order term, A = A^*'^ which is 
equivalent to assuming that A has no angle dependence. 
The equation is now expressed as 

A(p) = A / kdk^^^=C^'\klp\ (II) 
2p J ^k^ + A(fc)2 

where the Coulomb kernel is defined as 



27r (l-|-a;2-|-2a;cos6lfc)i/2 



(12) 



This integral equation can be solved numerically by it- 
eration, but it is instructive to discuss first an analytical 
solution that is available when a simplified version of the 
kernel is taken, given by 



C(")(a;) = 



(1 - .t) + ^61(0;- 1), 

X 



(13) 



which has the correct asymptotic behavior at large and 
small arguments. The Coulomb kernel and its approxi- 
mate form, as well as higher order kernels that are ne- 
glected, are plotted in Fig. [5] for comparison. The ad- 




FIG. 2. Plots of the first three circular harmonics of the 
Coulomb kernel C^'^\k/p), Eq. (fT2)) . The corresponding 
simplified kernel given in Eq. (|13p is shown for comparison 
(dashed line). 



vantage of this approximate kernel is that it allows one 
to re-express the integral equation as a differential one. 
Taking two derivatives of Eq. (jlip with respect to p, we 
obtain 



P 



9p2 



2p 



dA p pA 



0. 



(14) 



The boundary conditions for this equation are obtained 
from the first derivatives of Eq. (|TT|) at the endpoints of 
the integration region 



dA 

dp 



P 



FA 

dA 
dp 



= 0, 



0, 



(15) 
(16) 



where fi is an infrared cutoff (the inverse system size). 
Eq. (jl4p cannot be solved analytically, but one may lin- 
earize it to obtain the branching points in (3 where non- 
trivial solutions become possible^ii^, 



P 



,d^ 

dp"^ 



+ 2p 



dA p 



dp 



+ ^A = 0. 



(17) 



This equation is the well known Schroedinger equation 
with a l/r2 potentiaP^. Its solution is given by 



A = + A^p 



(18) 
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with i]± = 1/2(1 ± 7), 7 = y/1 — 2/3 and all quantities 
with a tilde are scaled with the cutoff p — p/A. Applying 
the infrared condition (Eq. [T5| we get 



(19) 



The ultraviolet condition (Eq. llSp is also homogeneous 
and so does not fix the scale of the gap A, something 
that can only be done with the non-linear equation. The 
ultraviolet condition gives is the set of branching points 
where non-vanishing solutions of the gap equation be- 
come possible 



a/2/3 — 1 log fl — 2 arctan 



1 



+ 2n7r, (20) 



which represents an infinite number of logarithmically 
spaced solutions /3c*^for 13 > 1/2. (When /i — ^ all so- 
lutions collapse to /3c = 1/2, although there is still one 
of lowest energy^,) . Note however that for (3 greater than 

I3c^\ a finite gap is always present and the linearized 
equation cannot be used. The rest of the solutions thus 
represent higher order instabilities that would take place 
if the first one is set artificially to zero. Coming back 
to the full integral equation, Eq. (ITT|) . we can solve it 
by iteration procedures, and obtain both the gap and its 
dependence on /3, which are plotted in Fig. [31 We ob- 
serve that, for fl = 10"^", a finite gap is generated for 
/3 > /3c*''' ~ 0.48, which is the signature of the excitonic 
transition. From now on we will simply call j3c = I3^c\ 
the physically relevant critical coupling. 
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FIG. 3. Dependence of the gap A(0) on /3 for /i = 10"^°. 
The gap becomes significant within numerical precision for 
/? > 0.48. Note that this computation is done with the full 
kernel C^"' and p'^^ need not be greater than 1/2. Inset: 
Momentum dependence of the gap A(p), for /3 = 0.61 (lower 
curve) and /? = 0.65 (upper curve), computed with Eq. (jlip . 



III. RESPONSE FUNCTIONS 

A complementary approach to study the excitonic 
transition is through response functions. These are de- 
fined as the expectation value of a mass operator Mip 



expanded to first order in an external mass perturbation. 
They can be considered as generalized susceptibilities in 
the mass channel, and their singularities represent the 
instabilities of the system to spontaneous mass genera- 
tion. Therefore they represent an alternative approach 
to study the excitonic transition. However, the impor- 
tance of the response functions goes beyond the char- 
acterization of the transition, as they represent physical 
observables with non-trivial interaction signatures also 
for /3 < /3c. This has also been emphasized in the con- 
text of the QED3 theory of the cuprates^i^^, where the 
behavior of the mass susceptibilities is similar to what is 
found in graphene. 

The response in the Af„ channel to a perturbation Mm 
is simply the correlator 



^ tr [M„G(p)r^(p,p + q)Gip + q)] , 

(21) 



which is shown diagrammatically in Fig. El^b) , and where 
r„i is the full mass vertex (note that Fm is a 4x4 matrix; 
the sublattice / valley index is omitted for clarity) and the 
factor of 2 accounts for spin. Again, this equation is exact 
but to solve it we need to approximate F in some way. To 
guide us in the choice of approximation, we realize that 
the full mass vertex satisfies the following property: if we 
include an external mass ninip'' Mnip in the Hamiltonian, 
then by construction it holds that 



dG{p) 



drrir. 



= G{p)r^{p,p)G{p)l 



(22) 



or equivalently in terms of the self-energy and writing the 
indices explicitly 



(M„)i 



dj:,,{p) 



dr 



= {^n)tj{p,p)\. 



(23) 



This identity can be proven diagrammatically by real- 
izing the derivative acts by cutting all possible fermion 
lines in the full self-energy and introducing a mass ver- 
tex Mn at every cut. If we assumed the rainbow sum- 
mation for the self-energy, it can be seen that the corre- 
sponding approximation for F is the ladder approxima- 
tion, depicted in Fig. Ufa). In this case F satisfies the 
self-consistent equation 



d^k G{k)Trn{k,q)G{k + q) 



\P- k\ 



(24) 

It is important to note that the propagators in this equa- 
tion are the full propagators in the rainbow approxima- 
tion. These are just Gq for /3 < /3c as shown in the pre- 
vious section, but acquire a mass for /3 > /3c and this has 
to be included to have a consistent computation. This 
simply represents the fact that after the phase transi- 
tion, the propagators have to be computed in the broken 
symmetry state. Note also that, as in what happens in 
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the gap equation, we anticipate by symmetry that the 
response function will satisfy — ^„„n, so we can 

drop the indices in this case too. 




<3> 



(b) 



FIG. 4. (a) Diagrammatic equation for the three point ver- 
tex (shaded triangle) in the ladder approximation. The cross 
denotes a mass vertex, (b) Response function diagram. 



We will now compute the response and vertex func- 
tions, and for the sake of clarity we will do so in two steps. 
To have some insight into the excitonic transition itself, 
it is simplest to discuss just Tl{qo = 0, 5 = 0). In addition 
the computation for /3 > /3c presents no major difficulty. 
Our main interest, however, is the full response function 
11(01,(7), which will compute next in what appears to be 
the experimentally relevant case of /3 < /3c. The exten- 
sion of n(w, q) to P > is also feasible but it is beyond 
the scope of this work. 



A. n(0, 0) and the excitonic transition 

The computation of 11(0, 0) is simplified by the identity 
in Eq. ()23|) . Since in the previous section we assumed 



= S„(M„)ij, Eq. (|23p implies that that r„(p,p) 
is also proportional to the corresponding mass matrix, 
(Xn)ij = {Mn)ijT. This allows us to perform the traces 
easily in both Eqs. (PT|) and ([M)) . Plugging in G as ob- 
tained in the previous section, and integrating in energy, 
the equations read 



n(g) 



(27r)2 (p2 + A2)3/2 



(2^)2 (fc2+A2)3/2 



(25) 

T{k,k). 
(26) 



As for the gap equation, we now assume that F has no 
angle dependence (i.e. we take only the first order of the 
expansion in circular harmonics) we find 



p 



(p2 + A2)3/2 

2^ 



dk 



(fc2 + A2)3/2 



r{p,p), (27) 
C^"\k/p)T{k,k). 



(28) 



A direct computation shows that Eq. ((28)) can also be 
obtained by adding an external mass m„ to the propaga- 
tor of Eq. (IH) and taking the derivative with respect to 
it, as Eq. mandates. 



We are now ready to make the equivalence between 
the gap equation and the response function approaches 
explicit. Consider Eq. ([SS)) for A = 0. If we write the 
equation for F in the form 



/ 



dkA{k,p)r{k,k) = 1, 



then 



A{k,p)^S{k-p)-fc^°\k/p). 

2p 



(29) 



(30) 



We can obtain F by inverting the operator A{k,p), i.e. 
F — A^^ . On the other hand, the linearized version of 
the gap equation has the form 



dkA{k,p)A{k) = 0, 



(31) 



in terms of the same operator. Therefore, whenever there 
is a non- vanishing solution of the linearized gap equation 
A(k,p) develops a zero eigenvalue, its inverse becomes 
singular, and F develops a divergence. Therefore the re- 
sponse function 11(0, 0), which is the just the integral of 

F, also develops divergences at the /3c''' whenever the lin- 
earized gap equation has a solution, i.e. at the critical 
points for the different instabilities. 

The computation of 11(0, 0) with A set to zero was 
carried out in Ref. [isl . where it was proven that an in- 
finite number of logarithmically spaced poles appear for 
/3 > /3c. This result is reproduced in Fig. [5] More- 
over, the analytical solution in terms of the model kernel 
showed that these poles are in fact given also by our Eq. 
(PO)) for the gap equation, as they should. The artificial 
constraint A = may therefore be used to locate the 
critical couplings for higher order solutions of the gap 
equation. Physically, however, as long as we cross /3c 
and the gap is generated, the system always stays in the 
lowest energy ground state and no further poles should 
be observed. This is indeed what is obtained from the 
numerical solution of Eq. (|28p if A is included as com- 
puted from Eq. pT|) . In Fig. [5] we show both cases for 
comparison. 



B. General computation of F and 11 

We now proceed to solve the general equations (PT|) 
and We will only consider /3 < /3c for simplicity, 

i.e., we set A = for the rest of the article. To solve this 
set of equations, it is convenient to decompose Tm in a 
basis of 4x4 matrices with well defined transformation 
properties under the SU(2) valley symmetry. Defining 
M = Tz(Tz, this basis may be taken as the four matrices 
Af ,X, a' which are scalars under this symmetry, and the 
matrices Mn,Tn, a^Tn, each of which transforms like a 
spin 1/2. With this choice we express F^ as 



F,„ = f,„M + f? T + fla' 



V" T 

^ nm-^ 11 



(32) 
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FIG. 5. Zero momentum response function 11(0,0) as a func- 
tion of /?. When the gap is not included in the electron prop- 
agators, the result presents logarithmically spaced poles (dot- 
ted black line). Asymptotes are shown (thin dotted hues) 
where poles occur. When self-energy diagrams are included 
in the fermion propagators, the generation of a mass prevents 
the appearance of poles after the transition (full red line). 
Note the first pole is approximately located at the critical 
coupling obtained from the gap equation. 



The equations are further simplified when is ex- 

pressed in terms of its longitudinal and transverse parts 



where q — q/ q. A similar relation applies for rjj,^. With 
the identities 



k ■ r„m — k ■ q r„„ k X q r„„. 



k X Tnm = k ■ q r„„ + kx qT 



nm ' 



(34) 
(35) 



substituting Eq. (15^ into Eq. (HI]), and performing the 
trace, we obtain 



n„m(g) 



d^p 8 



T 

nm 



(27r)3 D 

+pxg'(/i3r^„ + /i4r^„)], (36) 



where we have defined the denominator 

D{p,q) ^[pl-f + ze][(po + go)' - {P + q)^ + le], (37) 

and where all fij {p, q) (specified below) are even func- 
tions under the reversal of the relative angle 9p^^ — 
9p — 9q. Because of the decomposition in Eq. ([32]) , the 
scalar parts decouple completely and are not needed. We 
can then obtain equations for the relevant components of 
Tm by multiplying Eq. (1241) by the corresponding basis 



matrices and taking the trace. One then obtains 



(Fk 1 1 



(2^)2 D k 

+k X g(/i3r,^„ - /i4r°,j 

d^k 1 1 



[/nr. 



+ /i2r 

nm 



= ~-il3 



(2^)2 D Ip-^ k 

+ k X g (/23rL + /24r°„)) 



[/2ir„ 



/22r 



T 

nm 



(38) 



(39) 



r^j„ and r°„„ satisfy similar equations, but are not 
needed in what follows. We now perform a circular har- 
monic expansion 



p(np,ng) 



2tt 



2tt 



r(p,g), (40) 



and retain only the first order contribution. Terms con- 
taining k X q are odd and vanish. Thus, and F° 
completely decouple to first order. As anticipated, from 
the structure of Eqs. ([55]) and ([55]) it can be seen 

that in fact Ilnm = Snm^- With this simplification the 
relevant components of fij are 



fii 

fl2 
/22 

Defining 



-ko{ko + qa) + k{k + q), 

f y gO^g ; ^ 
J2l = i[— koq), 

2(9 X k)^ , /, N rn -^ 
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d9. 



fi'j 

D ' 



(41) 
(42) 

(43) 
(44) 



the self-consistent equations to first order in the circular 
harmonic expansion finally read 



r(o^n) = 1 



.(0.0) 



2p 
2p 



^(0)p(0,0) 



^12 ^ T 



(45) 



"22 ^ T /' 



(46) 



where the Coulomb kernel C'"' was defined in Eq. ^ 
The mixing Kernel if 12, as well as the higher order har- 
monics of the kernel Ku can be shown to be small and 
may be neglected. In this case, the final equations deter- 
mining the response function, spelling momenta explic- 
itly, are 

F(°.0)(p,g)^l + |- J dkC(°\k/p)K[f{k,q)r(°-"Hk,q), 

(47) 

U{q)^- f dpK['i\p,q)T("'"Hp,q). (48) 



7 



When the external q < qo, all Kij develop an imaginary 
part for {q^ ~ q)/2 < k < {qQ + q)/2. Note that when 

a; = g = 0, Kf^{k,q) — 1 and we recover Eq. ([28|l for 
A = 0. 



Analytic solution for the static vertex and 
response function for /3 < /3c 



In this section, we show how equations (|47l) and 
can be solved analytically in the static limit go = 0, if one 
assumes simplified versions for the kernels in the spirit 
of the previous section. This type of solution is also re- 
lated to the one employed in the Lorentz invariant case in 
QED3 in the computation of the propagator of fermion- 
antifermion composites^l (the analog of excitons in our 
case). The inclusion of an infrared cutoff in this approx- 
imation makes it excessively complicated, so we will set 
/i = for this section. 

The explicit expression for the static kernel in Eq. (|47)) 

is 



K[-\k/q) 



dOk 6™^*=/!: / kq q{k + q) 
'^'^ f + 2qk\k \k + q\ 



(49) 



Note we can write K[^{k,q) = K[^{k/q) when qo — 0. 
We will use the following simplified version 



(50) 



which is compared with the actual kernel and its higher 
order harmonics that are neglected in Fig. [6] 




FIG. 6. Plots of the first three circular harmonics of the kernel 
k[^i (k / q) , Eq. The corresponding simplified kernel ((SD} 

is shown for comparison (dashed line). 



We define F = F^"'"' for convenience. Its integral equa- 
tion in terms of the simplified kernels is 

F(p,g)=l + A r dk-C^'^\k/p)T{k,q) 
2p Jo 



q 



/ dkC^°\k/p)r{k,q), (51) 
J q 



which can be seen to reduce to the usual equation when 
q — )■ 0. The advantage of the simplified kernel k[^\ as 
we now show, is that we can separate F into its g > 
and q < parts 

F(fc,g) = r<{k,q)9{q - k) + r>{k,q)9{k - q), (52) 

and obtain two equations that can be solved separately 
and then matched. The equation for F^(p, g) is simply 
Eq. ([5T|) when p < q 



(3 f\ k (e{p-k) 6{k-p) 



r<(p,,)^i4r..if^ 

2 Jo 1 \ P 
+ f J dk^r>{k,q), 



while for p > q we have 

F>(p,<z) = l + f r dk-V<{k,q) 
2 Jo PI 

2 7q V P 



r<(fc,<7) 



(53) 



If we define 



ai{q) = I dk—T<{k,q), 
Jo 1 
.A ^ 

a2{q) = J dk-T>{k,q), 



T>{k,q). (54) 

(55) 
(56) 



we can write the two equations as 

r<{p,q) = l + ^a,{q) 

+^ r dk- + ^(fe-p) 

2 Jo 9 V P k 



T>{p,q) = l + ^lai(g) 
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dk 



Pq 
2p 

e{p-k) e{k-p) 



F<(fc,q), (57) 



F>(fc,g). (58) 



,q \ P « 

Taking two derivatives with respect to p we obtain 



P 



, a^F< 

dp^ 



2p- 



Pp 

dp ^ 2q 



F< = 0, 



p2^!£l _^2p— + ^F 
9p2 dp 2 



0. 



(59) 
(60) 



Thus we have obtained two independent equations for 
F< and F>. These equations are only coupled through 
the boundary conditions, which can be obtained from the 
first derivatives of equations ([57| and ([58|). For F< these 
are 



dp 



1 a2 

2 

^ dp 



0, 



= 0, 



(61) 
(62) 
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while for they are 



P- 



dp 



r> -1 



p- 



dp 



= 0, 



0. 



(63) 
(64) 



The differential equations have straightforward solu- 
tions. We note the equation for r> is again the same as 
the one obtained in the gap equation, and the one ob- 
tained in Ref. [H, with q playing the role of the infrared 
cutoff. The equation for is a Bessel-type equation. 
Their solutions are 



r> = A+p'^+ + A^p"^- 

\ -1/2 



(65) 



(66) 

Applying the boundary conditions we get the solutions 



(g"- -ai(l + 7]-)^) (l-t-77+)p"+ 



+ ('7+ ^7-), 
(67) 



(l + /3a2/2)Ji(v/Wg) fp 



-1/2 



(68) 



These solutions still depend on ai^2- Plugging them into 
Eqs. ([55)) and ([56]) we obtain a linear system of equations 
for a 1,2, whose solutions are 

(2//3 - 1) ((1 + 77+)2 - qf{l + r,.)^) + jqn+ 
((1 + 77+)2 - g7(l + ^_)2) /^(^) + /32/4(^~, _ 1) ' 

(69) 



ai 



0-2 = 4>{P)ai 



2 



with 



The response function can be finally obtained as 
Ii{q)^- r dp^T<ip,q) + - rdpT>ip,q). 

TT Jo 9 Jn 



(70) 



(71) 



(72) 



Evaluating the integral, plugging the values of ai^2 and 
in the limit g << A, we finally obtain 



1 



- [1 + {jvm - 1) 9" 



(73) 



TT (1 + ry^ 

This result reproduces the power law behavior of n(w = 
0, q) found numerically in Ref. [l^ and shows analytically 
that the exponent is in fact given by 770 = 7 = -^/l — /3//?c 
with (3c = 1/2. This is the excitonic transition again, in 
limit /i — >■ 0: when /3 = /3c the response function becomes 
singular. This analytical expression provides a simple, 
compact expression for observables that couple to the 
static response function in the ladder approximation. 



D. Dynamic response 

The kernels in the general frequency dependent re- 
sponse are too complicated for an analytic solution. 
Therefore, we now solve Eq. (I47p numerically by dis- 
cretizing the momentum fc on a logarithmic mesh and 
solving the corresponding matrix equation by Gaussian 
elimination. The integration of Eq. (|48|) is straight- 
forward. The results of this procedure are shown in 
Fig. [T] The self-energy is represented as the difference 
An = Il{qo,qo + 6q) - Tl{qo,qo) with Sq ^ qo ~ q for 
convenience. 
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<1 
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FIG. 7. Response function An (go, go -I- (Jg). a) An for \5q\ « 
qo and /3 — 0.36, real part (full line) and imaginary part 
(dashed line). Inset: the Kekule phonon displacements, b) 
Logarithmic plot of AH for /3 = 0.24 and Sq > (dotted line). 
The full lines are linear fits with 77 = 0.45 for Sq « qo and 
770 = 0.69 for 5q » qo- 



The main result of the inclusion of finite frequency is 
that the response is strongly modified at q close to qo, 
but remains essentially the static result once q » qo- 
Fig. [ZKa) displays the real and imaginary parts of An 
for \Sq\ « qo- We observe a cusp a.t 6q = in the 
real part, and a finite imaginary part for Sq < 0. Log 
plots of both sides of the real part and and the imaginary 
part reveal power laws as Sq ^ 0. A Kramers-Kronig 
analysis for \Sq\ << qo shows that this is only consistent 
if An cx (Sq)"^, i.e. the exponents are all the same^. Fig 
WO^) shows a log plot for 5q > where power law behavior 
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is evident for Sq << qg. We also observe that AH crosses 
over to a different power law for Sq >> qo, which we 
identify as the static result q^°^^, obtained analytically 
in the previous section. The inset of Fig[7l^b) shows that rj 
is /3-dependent, and that it tends to the non- interacting 
result in Eq. (|80|) as /3 — ?> 0. In summary, the main 
features of the response function are a cusp (g — go)'' 
around go, and a crossover to the static power g''° for 
q » go- 



IV. EXPERIMENTAL SIGNATURES OF 
POWER LAW BEHAVIOR 

In the previous section we found that the mass sus- 
ceptibilities have characteristic power law behavior with 
/3 dependent exponents. We may now ask what are the 
experimental consequences of this. In general, these cor- 
relations can be observed in a linear response-type exper- 
iment, with a suitable probe that couples to electrons in 
the form of a mass. This is however difficult, as the usual 
experiments rather couple to the electron charge or cur- 
rent. To find a probe that couples to the masses we need 
to refer to their microscopic origins: A Kekule distortion 
for Mi,M2, or a sublattice antisymmetric potential for 

Af3. 

A first proposal to measure the II33 correlator was put 
forward in Ref. [l^ which involved placing a Coulomb 
impurity asymmetrically with respect to the sublattices, 
and measuring the sublattice charge difference with a 
STM tip. This measurement is difficult to perform, as 
it requires one to resolve the lattice structure in detail. 
A different possibility that we now discuss is to probe 
particular phonons that couple to electrons with a mass 
vertex. The self-energy of this type of phonons is pre- 
cisely given by the mass susceptibility, which then be- 
comes observable through the dispersion and lifetime of 
the phonon. These can be measured with current exper- 
imental techniques discussed below. 

The phonon spectrum of the honeycomb lattice con- 
sists of six phonon branches, four in-plane and two out- 
of-plane. Each of these phonons may couple to elec- 
trons near either Dirac point if it has momentum close 
to zero (a F point or zone center phonon), which scat- 
ters electrons within each valley, or if it has momen- 
tum close to K or K' points (a zone boundary phonon), 
in which case it produces intervalley scattering. The 
strength of the electron- phonon coupling (EPC), how- 
ever, depends on how the particular displacement pattern 
of that phonon modifies the hopping integrals between 
atoms. Two modes have displacements that produce a 
significant EPC, and both of them are in-plane phonons. 
The first of these is the phonon branch of highest energy 
at the F point, the E2 phonon. The second is the Ai 
branch at the K and K' points (also the highest branch). 
This is a lattice distortion with a supercell of six atoms, 
whose displacement pattern is obtained by taking linear 
combinations of the displacements at K and K' , and is 



shown in the inset of Fig. |8l These two combinations 
couple to electrons exactly in the same way as the two 
components of the Kekule distortion, i.e. they couple 
with the mass matrices Mi and M2 

He-ph,K = Pk J dPrip\MiUKi + MaMAra)^', (74) 

with Fk — idt/da. For this reason this phonon is also 
known as the Kekule phonon^^ . This phonon is therefore 
suitable to measure the IIn — 1122 correlators. 

It has also been noted that in the presence of broken 
z ^ —z symmetry, induced for example by a substrate or 
a perpendicular electric field, there is another mechanism 
that produces an EPC with out-of-plane phonons^^i^. 
This is simply that atoms displaced to different positions 
see a different potential, and is not related to changes 
in hopping integrals. In particular, for the ZO (out-of- 
plane, optical) phonon at the F point, the sublattices A 
and B move in opposite directions, and there is a linear 
coupling of the form 



ph.ZO 



Fzo / d^r%p'^ M^uzoi', 



(75) 



with Fzo Ez- Therefore, in the presence of an electric 
field, the ZO phonon can be used to measure 1133. The 
Hamiltonian of either phonon may be expressed as 



H 



E 



<fq 
(2^ 



(76) 



with creation and destruction operators defined by 



Ac 



ifq 

AuJiM J (27r)2 



{h,,e^^'+by-^^~), {11) 



where i = Kl, K2, ZO, ujk « 0.17 eV, ujzo ~ 0.1 eV, A^ 
is the unit cell area, and the dispersion of the phonons 
is neglected to a first approximation. A dimensionless 
EPC can be defined as A^^ = F^Ac/{2Muj,vl). In the 
case of the Ai phonon it is estimated to be in the range 
Xk ~ 0.03 - 0.1^^. We will now discuss the only the 
Ai phonon, as the ZO case in an electric field has the 
same behavior. 

The phonon propagator can be obtained in terms of 
the self-energy as 



Gph{uj,q) = 



2lok 



(78) 



and it follows directly from the form of the electron- 
phonon vertex Eq. (I74p that the phonon self-energy is 
directly related to the mass susceptibility Eq. (|211) via 
S = The dispersion relation for the phonon can 

be obtained by solving for the pole in Eq. ([75)) for small 
\k , so that the dispersion relation is corrected to 



a;(g) w Wif ^n(wif, g) 



(79) 
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For /3 = this isSi 

co{q)^UK + ^{vW-^]<y^', (80) 

which has a square root singularity at qk for q > qx- For 
q < qK the self-energy is purely imaginary, and a finite 
lifetime is obtained. The Kohn anomaly is conventionally 
associated with a linear cusp in the dispersion, which is 
obtained only asymptotically for q » qK', the full dy- 
namical self-energy should be used in general. Note that 
qx is approximately 2% of the T — K distance in the 
Brillouin zone. The necessity of employing the dynami- 
cal self-energy has been emphasized before^^— , in par- 
ticular in the doped case where the static approximation 
produces poor agreement with experiments^. Note also 
that a different Dirac fermion induced Kohn anomaly has 
been recently observed in the surface of topological insu- 
lator BisSe^i^. 

In the presence of electron-electron interactions, E is 
modified as described in the previous section and in par- 
ticular acquires power law behavior. In Fig. |8] we plot 
the phonon dispersion relation obtained from it for differ- 
ent values of /3. This is given in terms of the self-energy 
evaluated at the phonon frequency luk ■ To ease the com- 
parison at different values of /3, we also represent the 
difference 

Aa;(g) = uj{q) - uj{qK) = ^ (n(wK, q) - n(wK, qx)) , 

(81) 

where we have recovered physical units with hvp — 6.5 
eVA. The values of the parameters used are Xk =0.1 
and A — 1.7eV. The dispersion follows the static power 
law q^oiP) for q >> q^:, and the cusp turns into q'^^^^ as 
discussed above. The modification of the Kohn anomaly 
due to interactions is rather dramatic and should be ob- 
servable. 




0.00 0.05 0.10 0.15 0.20 0.25 

q(A-l) 

FIG. 8. Ai phonon dispersion relation ALu{q) measured from 
the K point for /3 — 0,0.1,0.2,0.3, with higher curves cor- 
responding to higher values of j3. Note that a;(gx), which 
depends on /3, has been subtracted from each curve for an 
easier comparison. Inset: the Kekule phonon displacements. 



From the experimental point of view, there are sev- 
eral techniques available for the measurement of the Ai 
phonon dispersion, and each one has its own potential dif- 
ficulties. One method is to employ Electron Energy Loss 
Spectroscopy (EELS), which has already been used to 
map K-point phonon dispersions of graphene on different 
substrates where it behaves as quasi-freestanding^^^, 
such as Ft. (The absence of hybridization with the sub- 
strate is important as it would strongly change the elec- 
tron band structure and the Kohn anomaly^) . The prob- 
lem with metallic substrates is that screening will turn 
the long range Coulomb interaction between graphene 
electrons into a short-ranged one, or may even suppress 
it completely. Thus a power law anomaly is not expected, 
which is consistent with the linear one observed in Ref. 
H . An insulating substrate would be required to observe 
this effect. 

A more indirect experiment (with insulating substrate) 
is to track the dependence of the 2D Raman peak with 
incoming laser energy. This method has been used^ to 
measure the dispersion of the Ai phonon, but both the 
amount of data it yields and the range of momenta it cov- 
ers is limited and not very close to the K point. Finally, 
X-rays are a usual tool to measure phonon dispersions 
in 3D crystals, and while it is probably challenging to 
obtain enough intensity from a single sheet of graphene, 
experiments in graphite^^ii^ might be used to deduce the 
phonon dispersion. This approach is not straightforward 
because the electronic structure of graphite is different 
from graphene, and this must be taken into account. 
Nevertheless, it is encouraging to observe that precision 
measurements show an Ai phonon dispersion that is not 
linear^®. 



V. CONCLUSIONS AND DISCUSSION 

Understanding the role of interactions in graphene is a 
challenging problem whose solution is far from complete. 
One of the main purposes of this work is to bring greater 
attention to a new potential way of observing interac- 
tion effects: the measurement of the mass susceptibil- 
ity, which, as we have shown, displays power law behav- 
ior with /^-dependent exponents. Apart for unveiling a 
novel signature of interactions, the identification of these 
power laws could represent an alternative measurement 
of graphene's fine structure constant /3. 

Moreover, the observation of the static power law 
would also help in understanding the problem of the ex- 
citonic transition in graphene. With the experimental 
evidence gathered so far, there seems to be no indica- 
tion of the presence of an interaction induced gap, even 
in high quality samples^i. Theoretically, however, it ap- 
pears that suspended graphene should be, if not in the 
gapped phase, at least close to the transition. Since, as 
we have shown, the static exponent of the mass-mass 
response goes to zero at the critical /?c, a measurement 
of the exponent would indicate how close we are to the 
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achievement of the long-sought gap. 

Note also that other potential measurements of these 
power laws could come from spin related experiments, 
which we have not discussed. The inclusion of spin in 
the picture is straightforward: in this case, the symme- 
try of the Hamiltonian enlarged to SU(4), and the ex- 
citonic masses generalize to singlet and triplet versions. 
Because of the same symmetry argument relating CDW 
and Kekule masses, any of these correlators would have 
the same power law behavior, so the observation of these 
particular spin susceptibilities may represent an alterna- 
tive route to the experiments we have discussed in this 
work. 

A final comment concerns the robustness of our result 
to more refined approximations schemes than the ladder 
summation. While other sets of diagrams may modify 
our quantitative predictions, it is very unlikely that the 
non-analytic behavior can be removed in this way. One 
may consider, for example, the inclusion of self-energy 
terms for the electron propagator—, which may produce a 



slow logarithmic dependence of the exponent. This could 
be taken into account in a similar way as was done for the 
gap equationi^. The effect of static RPA screening of the 
Coulomb interaction will in general reduce the effective 
value of f3 but will not change the power law itself. Fi- 
nally we also note that the 1/N approximation does gives 
power law behavior for the Kekule mass correlato r^"'^^ 
(and thus the self-energy) as well. 

In summary, we believe that the observation of power 
law correlations originating in the mass susceptibility is 
potentially feasible and would be an important contribu- 
tion to our understanding of the problem of interactions 
in graphene. 
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